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ABSTRACT 

In this paper an ocean scatter propagation 
model, developed for aircraft-to-satellite 
(aeronautical) applications, is described. The 
purpose of the propagation model is to 
characterize the behaviour of sea reflected 
multipath as a function of physical propagation 
path parameters. An accurate validation against 
the theoretical far-field solution for a perfectly 
conducting sinusoidal surface is provided. 
Simulation results for typical L-band aeronautical 
applications with low complexity antennas are 
presented. 

INTRODUCTION 

For L-band oceanic mobile satellite 
communications, sea reflected multipath is the 
most significant propagation component, 
especially when relatively low gain, non-directive 
antennas are employed. A propagation model is 
required to characterize the behaviour of sea 
reflected multipath as a function of the elevation 
angle to the satellite, aircraft altitude, sea state 
(surface roughness), aircraft antenna 
characteristics, aircraft velocity, and signal 
polarization. 

Evaluating the electromagnetic field scattered 
by a rough surface is an extremely difficult 
problem. Most theoretical approaches are based 
on the Helmholtz surface integral, with the 
electric field on the surface assumed to match that 
present on a tangent plane at that point [1]. In 
general, complicated integrals that are difficult to 
solve are obtained. 


In many theoretical treatments, attention is 
restricted to perfectly conducting surfaces in 
order to simplify the scattering solution. In the 
perfectly conducting case, the surface reflection 
coefficient has unity magnitude, independent of 
grazing angle. This is not appropriate for the sea 
surface, as is evident from the behaviour of the 
reflection coefficients at L-band shown in Figure 
1. In most theoretical treatments, shadowing and 
blockage by the surface is ignored for 
mathematical tractability. Unfortunately, this 
makes low elevation angle scatter predictions 
suspect. The propagation model described in this 
paper incorporates the variability of the surface 
reflection coefficient as well as shadowing and 
blockage effects. 



Figure 1: Sea surface reflection^ coef- 
ficients at L— band ( A c — 0.2 m ). 
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DESCRIPTION OF THE PHYSICAL 
PROPAGATION MODEL 

A simplified view of the physical model that 
was implemented as a computer simulation is 
provided in Figure 2. Along the azimuthal 
direction to the satellite, a rough wave profile 
z w (x) is superimposed on a curved earth surface, 

x 2 

z s (x) = z w (x) - ^ (1) 

where ae is the radius of the earth. 

It is assumed the waveheight does not vary 
perpendicular to the azimuthal direction. The 
wave profile is described by a sum of random 
phase sinusoidal components of sufficient 
number to give reasonable agreement with 
tabulated sea state data. A summary of the wave 
profile model and appropriate sea state indices is 
provided in Appendix A. 

The initial task in the determination of the 
scatter signal received at the aircraft antenna is the 
identification of the surface regions making the 
most significant contributions. These are the 
regions about the "specular" points, which 
geometrically reflect rays to the receiver (see 
Figure 2). Locating specular points involves 


finding the values of x where 

0 (x) = B (x) (2) 

with 

0 (x) = tan 1 Zs(x) I (3) 

[ A - A 0 

B (x) = E - 2 £(x) (4) 

£ (x) = tan' 1 (5) 


. . E + B(x) , c , 

a(x) = E - e(x) = ~ * (6) 

Here, £(x) denotes the incline of the tangent plane 
to the surface, B(x) is the angle of the reflected 

ray, a(x) is the local grazing angle, the co- 
ordinates of the antenna are (xo,zo) and E is the 
elevation angle. The specular point search 
procedure is extremely computationally 
demanding for aeronautical applications (i.e. the 
active scattering region extends to the horizon 
range of the aircraft). An attractive feature of the 
simulation model is the exclusion of scattering 
facets where the specular point is blocked by the 
wave profile. 
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In the propagation simulation, unmodulated 
transmissions are considered and the H-pol and 
V-pol components are treated separately. The 
scatter signal from each contributing surface facet 
is accurately determined. This is accomplished by 
incorporating the variability of the surface 
reflection coefficient, aircraft antenna gains, as 
well as the divergence and phase shift resulting 
from scattering from curved surface facets. That 
is, 

S H (i) = r H ( aO e-j^i Gh(B A - Bi) Ai ( 7 ) 

S V (i) = T V ( «;) e'i 0 i G v (B a - Bi) Ai (8) 

where Th, Tv and Gh, Gy denote the surface 
reflection coefficients and mobile antenna gains 
for H-pol and V-pol, Ba is the aircraft banking 
angle (positive towards the surface), 0 i is the 
scatter path phase shift (relative to the direct 
path), Aj is the signal attenuation and phase shift 

factor, and B, and cij are the reflection and 
grazing angles, respectively. 

The scatter path phase shift can be expressed 
as 

0 _ 2tc rj (1 - cos(E-t-Bj)) 

Xc 

where rj is the distance between the antenna and 
the scattering facet and is the carrier 
wavelength. 

An appropriate signal attenuation and phase 
shift factor for the aeronautical application, 
derived in Appendix C of [2], is 


^1 + pj 1/2 for convex facets 
j - lV 1/2 for concave facets 


Here, fj is the focal length of the curved facet and 
r c is the surface radius of curvature. For 
aeronautical applications, the magnitude of Ai is 
equivalent to the standard ray-optics divergence 
factor. 

The overall scatter signals for each 
polarization are given by the superposition of the 
individual components from each contributing 
surface facet. These signals can be subsequently 
combined to account for circular and elliptical 
polarizations. To give some appreciation for the 
complexity of the aeronautical propagation 
simulation, there are in excess of 10,000 
contributing facets for an aircraft altitude of 9.1 
km, an elevation angle of 10 °, and rough surface 
conditions. 


AN IMPORTANT VALIDATION OF 
THE SIMULATION SOLUTION 

For perfectly conducting, sinusoidal surfaces, 
the far-field Helmholtz integral solution is 
provided in Section 4.3 of [ 1 ]. The received 
electric field can be expressed as a scaled version 
of the field that would have been present in the 
absence of surface roughness. The scale factor 

p is a function of the angle of incidence of the 
incoming radiation and the scattering direction 
under consideration. Restricting attention to the 
specular direction, the theoretical scale factor can 
be expressed as 


Plheory = J< 


where 0 <j is the phase difference between the 
crest and trough specular paths, and J 0 ( ) is the 
Bessel function of the first kind of order zero. 

The scale factor obtained with the simulation 
model (p.61 of [ 2 ]) is 


where 


fi = y sin(oq) 


(1 -t- Zs(Xj) 2 ) 3 / 2 
Z^(Xi) 


Psim — 


cos (^2 ' V) G 4 ) 


The excellent agreement between the theoretical 
and simulation solutions for rough surface 

conditions (0d > jjpjis evident in Figure 3. 
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SIMULATION RESULTS 

Some L-band (f c =1.5 GHz) aeronautical 
simulation results for a low complexity crossed- 
slot antenna are presented in this section. The 
antenna pattern considered is presented in Figure 
4. For rough surface conditions, it has been 
confirmed that the multipath process can be 



<Pd / 71 


Figure 3 : Simulation and theoretical 
reflection factors for perfectly 
conducting sinusoidal surfaces. 


0 = 90° 



Figure 4 : Crossed-slot antenna pattern 


accurately modelled as a Rayleigh fading process 
[2], Given the large number of contributing 
surface facets, this is not a surprising result. 
With the addition of a direct path component, a 
Rician channel model is appropriate for 
aeronautical applications. The average 
interference level and behaviour of the second 
order statistics of the multipath process are 
required to characterize the fading aeronautical 
channel . 

Table 1 summarizes results in terms of signal 
to interference ratio (S/i) as a function of elevation 
angle for the crossed-slot and reference omni- 
antenna patterns. The aircraft altitude was 9. 1 km 
and the surface conditions were rough (sea state 
index SS4(ii)). The superiority of circular 
polarization (C-pol) is evident in Table 1. This is 
a consequence of the polarization sense reversal 
that a surface imparts on a scattered signal, which 
is discriminated against by the receiving antenna. 

For a circularly polarized crossed-slot antenna, 
it is evident that S/i is not very sensitive to 
elevation angle, and the values are excellent, 
exceeding 20dB. At low elevation angles, the 
crossed-slot antenna discriminates against H-pol 
multipath (see Figure 4), while the surface 
discriminates against V-pol multipath (see Figure 
1). With S/i values this high, ocean scatter will 


Elevation 

Angle 

S/I IdBl For 1 

H-nol 

V-pol 

C-pol 

31° 

24.6 

(1.2) 

16.9 

(3.9) 

21.6 

(17.9) 

20“ 

14.9 

(1.0) 

15.7 

(5.9) 

20.4 

(13.6) 

15“ 

11.9 

(1.6) 

16.5 

(8.0) 

21.3 

(12.4) 

10“ 

10.4 

(2.7) 

18.8 

(11.3) 

22.9 

(11.7) 

5“ 

11.2 

(5,2) 

21.0 

(16.2) 

22.9 

(12.5) 


(Values in parenthesis for an omni-antenna) 

Table 1: S/i as a function of elevation angle 
for an aircraft altitude of 9. 1 km 
and SS4(ii) surface conditions. 
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not have a severe impact on data communications 
integrity. Even for an omni-antenna, the S/j 
values are higher than 10 dB for circular 
polarization. 

The simulation results presented in Table 2 
demonstrate that S/I is not very sensitive to sea 
state and antenna height. The largest difference 
between the SS3 and SS5(i) results is only 0.7 
dB, while the biggest discrepancy between the 
results for antenna heights of 9.1 km and 5 km is 
only 0.4 dB. A slight trend of increasing scatter 
power as the surface gets smoother (decreasing 
sea state) and as the aircraft altitude is lowered is 
noticed. 

The power spectrum of the multipath process, 
commonly referred to as the Doppler spectrum, 
conveys information about the second order 
statistics of the fading process. An example 
Doppler spectrum for an elevation angle of 10' is 
presented in Figure 5. The most striking feature 
is the asymmetry. This can be theoretically 
justified (Appendix E of [2]), although the 
explanation is too lengthy to be presented here. 
The following relationship between Doppler 
spread Dio (Hz), horizontal airspeed v x (m/s), 
and elevation angle E, was derived from 
simulation results presented in Section 3.2 of [2]: 


( 0.23 v x , E=10' 
D 10 = / 0.37 v x , E=15° 
1 1.15 v x , E=31° 


(15) 



(Values in parenthesis for an omni-antenna pattern) 


Here, Djo is the width of the spectral region 
where the relative Doppler response is above -10 
dB. For an aircraft speed of 600 mph, the 
Doppler spreads are around 60 Hz at low 
elevation angles (E=10°) and around 300 Hz for 
intermediate elevation angles (E=31°). 

Figure 6 compares simulated H-pol and V- 
pol scatter power levels for an omni-reference 



Normalized Frequency (f/f 3 ) 

Figure 5 : Doppler spectrum for E - 1C 
and a C — pci omni-antenna 



Table 2: Effect of sea state and aircraft altitude 
on S/i for E=10' 


Figure 6 : Measured and simulated scat- 
ter power levels for an omni-antenna. 
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antenna against measurement results from [3] and 
a simple approximation given by the magnitude 
squared value of the product of surface reflection 
coefficient and the divergence factor of the earth 
[3]. The agreement for elevation angles above 
10* is quite good. 


A nice feature of this surface model is that 
derivatives, which are required in the geometrical 
search for specular points and radius of 
curvature computations, are easily determined. 
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APPENDIX A 
SURFACE MODEL 


The wave profile used in this study is 
described by a sum of random phase sinusoidal 
components, with a white bandpass spatial 
spectrum. 


z w (x) = 


V2o z N 

VNt k--N 


COS (271 Uk x + 0k) 


where Nt=2n+ 1 is the number of sinusoidal 
components, 0k is a uniformity distributed 

ky 

random phase, Uk=u 0 (l+j^), y= 1.344, and 

u 0 = (1.0726 A) -1 . The dependence of the rms 

waveheight Oz and average sea wavelength A on 
sea state is provided in Table A.l for a few 
representative cases. 


For the simulation runs, Nx was chosen to be 
21. With the above parameter values, the wave 
profile generated by this model is approximately 

Gaussian with a correlation distance of A/4.44. 
According to [4], this is a realistic model for a 1- 
dimensional cross-section of the ocean surface. 
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Sea 

Slate 

Index 

RMS 

Waveheight 

(Meters) 

Average Sea 
Wavelength 
(Meters) 

3 

0.35 

22.3 

4(i) 

0.525 

30.75 

4(ii) 

0.61 

34.3 

5(0 

0.76 

41.9 

5(H) 

0.92 

48.8 


Table A.l: Representative sea state table 
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